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Abstract. 

I summarize the status of three-neutrino oscillations that follow from combining the relevant world's data. The discussion 

includes the small parameters a = A#ig OL /Am^ TM and sin 2 613, which characterize the strength of CP violation in neutrino 

oscillations, the impact of oscillation data on the prospects for probing the absolute scale of neutrino mass in j3j3o v and the 

robustness of the neutrino oscillation interpretation itself in the presence of non-standard physics. I also comment on the 

theoretical origin of neutrino mass, mentioning recent attemps to explain current oscillation data. 
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INTRODUCTION 

The discovery of neutrino oscillations has marked a turn- 
ing point in our understanding of nature and has brought 
neutrino physics to the center of attention of the particle, 
nuclear and astrophysics communities. Here I summarize 
the determination of neutrino mass and mixing parame- 
ters in neutrino oscillation studies following Ref. [1] to 
which the reader is referred for details on data analysis 
and experimental references. For future neutrino oscil- 
lation projects see Ref. [2]. The structure of the three- 
flavour lepton mixing matrix in various gauge theories 
of neutrino mass was given in [3]. Current neutrino os- 
cillation data are well described by its simplest unitary 
form, with is no sensitivity to CP violation. The effect of 
Dirac CP phases in oscillations and Majorana phases in 
Pq v constitute the main challenge for the future. The in- 
terpretation of the data requires good calculations of so- 
lar and atmospheric neutrino fluxes [4, 5], neutrino cross 
sections and experimental response functions, as well the 
inclusion of matter effects [6, 7] in the Sun and the Earth. 

SOLAR AND KAMLAND DATA 

The solar neutrino data includes the rates of the chlorine 
experiment (2.56 ± 0.16 ± 0.16 SNU), the results of the 
gallium experiments SAGE (66.9 tH ±H SNU) and 



GALLEX/GNO (69.3 ±4.1 ± 3.6 SNU), as well as the 
1496-day Super-K data (44 bins: 8 energy bins, 6 of 
which are further divided into 7 zenith angle bins). The 
SNO data include the data from the salt phase in the form 
of the neutral current (NC), charged current (CC) and 
elastic scattering (ES) fluxes, the 2002 spectral day/night 
data (17 energy bins for each day and night period) and 
the 391-day data. The analysis includes not only the 
statistical errors, but also systematic uncertainties such 
as those of the eight solar neutrino fluxes. 

KamLAND detects reactor anti-neutrinos at the 
Kamiokande site by the process v e + p — ► e + + n, where 
the delayed coincidence of the prompt energy from the 
positron and a characteristic gamma from the neutron 
capture allows an efficient reduction of backgrounds. 
Most of the incident v e 's come from nuclear plants at 
distances of 80 — 350 km from the detector, far enough to 
probe large mixing angle (LMA) oscillations. To avoid 
large uncertainties associated with the geo-neutrino flux 
an energy cut at 2.6 MeV prompt energy is applied for 
the oscillation analysis. 

The first KamLAND data correspond to a 162 ton- 
year exposure gave 54 anti-neutrino events in the final 
sample, after all cuts, while 86.8 ±5.6 events are pre- 
dicted for no oscillations with 0.95 ± 0.99 background 
events, consistent with the no-disappearance hypothesis 
at less than 0.05% probability. This gave the first evi- 



dence for the disappearance of reactor neutrinos before 
reaching the detector, and thus the first terrestrial con- 
firmation of oscillations with Amg OL . Additional Kam- 
LAND data with a somewhat larger fiducial volume of 
the detector were presented at Neutrino 2004, corre- 
sponding to an 766.3 ton-year exposure. In total 258 
events have been observed, versus 356.2 ±23.7 reactor 
neutrino events expected in the case of no disappearance 
and 7.5 ±1.3 background events. This leads to a confi- 
dence level of 99.995% for v e disappearance. Moreover 
evidence for spectral distortion consistent with oscilla- 
tions is obtained. 

A very convenient way to bin the latest KamLAND 
data is in terms of 1 /E pl , rather than the traditional bins 
of equal size in £p r . Various systematic errors associated 
to the neutrino fluxes, backgrounds, reactor fuel compo- 
sition and individual reactor powers, small matter effects, 
and improved v e flux parameterization are included [1]. 
This singles out the LMA solution from the previous 
"zoo" of alternatives [8]. The stronger evidence for spec- 
tral distortion in these data also leads to improved Am^ 
determination, substantially reducing the allowed region 
of oscillation parameters. From this point of view Kam- 
LAND has played a key role in the resolution of the solar 
neutrino problem. Assuming CPT invariance one can di- 
rectly compare the information obtained from solar neu- 
trino experiments with the KamLAND reactor results. 

ATMOSPHERIC AND K2K DATA 

The first evidence for neutrino oscillations was the zenith 
angle dependence of the jti-like atmospheric neutrino 
data from the Super-K experiment in 1998, an effect also 
seen in other atmospheric neutrino experiments. How- 
ever, though appealing, the original oscillation interpre- 
tation was certainly not unique [9]. Today, thanks to the 
accumulation of upgoing muon data, and the observation 
of the dip in the L/E distribution of the atmospheric 
survival probability, the signature for atmospheric neu- 
trino oscillations has become clear. The data include 



Super-K charged-current atmospheric neutrino events, 
with the e-like and jU-like data samples of sub- and multi- 
GeV contained events grouped into 10 zenith-angle bins, 
with 5 angular bins of stopping muons and 10 through- 
going bins of up-going muons. We do not use v T appear- 
ance, multi-ring jj, and neutral-current events, since an 
efficient Monte-Carlo simulation of these data would re- 
quire further details of the Super-K experiment, in par- 
ticular of the way the neutral-current signal is extracted 
from the data. We employ the latest three-dimensional 
atmospheric neutrino fluxes given in [5]. 

disappearance over a long-baseline probing the 
same Am 2 region relevant for atmospheric neutrinos is 
now available from the KEK to Kamioka (K2K) neutrino 
oscillation experiment. Neutrinos produced by a 12 GeV 
proton beam from the KEK proton synchrotron consist of 
98% muon neutrinos with a mean energy of 1 .3 GeV. The 
beam is controlled by a near detector 300 m away from 
the proton target. Comparing these near detector data 
with the content of the beam observed by the Super- 
K detector at a distance of 250 km gives information on 
neutrino oscillations. 

The data K2K-I sample (4.8 x 10 19 protons on target) 
gave 56 events in Super-K, whereas 80.1 4 were ex- 
pected for no oscillations. The K2K-II data correspond 
to 4.1 x 10 19 protons on target, comparable to the K2K- 
I sample. Altogether they give 108 events in Super-K, 
to be compared with 150.9^}q q expected for no oscilla- 
tions. Out of the 108 events 56 are so-called single-ring 
muon events. This data sample contains mainly muon 
events from the quasi-elastic scattering v M + p — > fl + n, 
and the reconstructed energy is closely related to the true 
neutrino energy. The K2K collaboration finds that the ob- 
served spectrum is consistent with the one expected for 
no oscillation only at a probability of 0. 1 1 %, whereas the 
best fit oscillation hypothesis spectrum has a probability 
of 52%. 

One finds that the neutrino mass-squared difference 
inferred from the disappearance in K2K agrees with 
atmospheric neutrino results, providing the first confir- 



mation of oscillations with Aw 2 TM with accelerator neu- 
trinos. Unfortunately in the current data sample K2K 
gives a rather weak constraint on the mixing angle, due 
to low statistics. However, although the determination of 
sin 2 ATM is completely dominated by atmospheric data, 
K2K data already start constraining the allowed Am 2 TM 
region [1]. In particular, there is a constraint on Am 2 TM 
from below, which is important for future long-baseline 
experiments, since these are drastically affected if Am^ TM 
lies in the lower part of the 3(7 range indicated by current 
atmospheric data alone. 



THREE-NEUTRINO OSCILLATIONS 

The first systematic study of the effective lepton mixing 
matrix in gauge theories of massive neutrinos was given 
in [3]. For some models this matrix can be taken as 
approximately unitary. For three neutrinos, this gives 



(1) 



where each factor is effectively 2x2 and contains an 
angle and a CP phase. Two of the three angles are in- 
volved in solar and atmospheric oscillations, so we set 
012 = ©sol and 023 = 9 a™- The last angle in the three- 
neutrino leptonic mixing matrix is 0i3, 
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for which only an upper bound currently exists. All three 
phases are physical [10], one corresponds to the one 
present in the quark sector (Dirac-phase) and affects neu- 
trino oscillations, while the other two are associated to 
the Majorana nature of neutrinos and show up in neutri- 
noless double beta decay and other lepton-number vio- 
lating processes, but not in conventional neutrino oscil- 
lations [10, 11]. 

Current neutrino oscillation experiments are insensi- 
tive to CP violation, thus we neglect all phases. In this 
approximation three-neutrino oscillations depend on the 
three mixing parameters sin 2 0i2, sin 2 023, sin 2 0^ and on 
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FIGURE 1. Three-neutrino regions allowed by the world's 
oscillation data at 90%, 95%, 99%, and 3ct C.L. for 2 d.o.f.In 
top panels A^ 2 is minimized wrt undisplayed parameters. 

the two mass-squared differences Am 2 OL = Am 2 j = m 2 — 
;«i and Am 2 TM = Am 2 t = m 2 — m 2 characterizing so- 
lar and atmospheric neutrinos. The hierarchy Am 2 OL <C 
Am 2 TM implies that one can set, to a good approximation, 
Am 2 OL = in the analysis of atmospheric and K2K data, 
and Am 2 TM to infinity in the analysis of solar and Kam- 
LAND data. Apart from the data already mentioned, the 
global oscillation analysis also includes the constraints 
from the CHOOZ and Palo Verde reactor experiments. 

The results of the global three-neutrino analy- 
sis are summarized in Fig. 1 and in Tab. 1, taken 
from Ref. [1]. In the upper panels of the figure 
the Ax 2 is shown as a function of the parameters 
sin 2 0i 2 , sin 2 023 , sin 2 0i 3 , Am\ i , Am 2 j , minimized with 
respect to the undisplayed parameters. The lower panels 
show two-dimensional projections of the allowed re- 
gions in the five-dimensional parameter space. The best 
fit values and the allowed 3(7 ranges of the oscillation 
parameters from the global data are summarized in 
Tab. 1. This table gives the current status of neutrino 
oscillation parameters. 

As it has long been noted, in a three-neutrino scheme 
CP violation disappears when two neutrinos become de- 
generate [3] or when one angle vanihes, such as 0i3 [12]. 
All genuine three-flavour effects involve the mass hierar- 
chy parameter a = Am 2 OL /Am 2 TM and the mixing angle 

013. 



TABLE 1. Current oscillation parameters. 



parameter 



best fit 



3o~ range 



Am^j [10" 


- 5 eV 2 ] 


7.9 


7.1-8.9 


Amfj [10~ 


- 3 eV 2 ] 


2.2 


1.4-3.3 


sin 2 On 




0.31 


0.24-0.40 


sin 2 023 




0.50 


0.34-0.68 


sin 2 0i3 




0.000 


< 0.047 





FIGURE 2. Determination of a = Am 2 OL /Am 2 TM and bound 
on sin 2 $13 from current neutrino oscillation data. 

The left panel in Fig. 2 gives the parameter a, namely 
the ratio of solar over atmospheric splittings, as deter- 
mined from the global % 2 analysis of [1]. 

The right panel in Fig. 2 gives A.% 2 as a function of 
sin 2 0i3 for different data samples. One finds that the 
KamLAND-2004 data have a surprisingly strong im- 
pact on this bound. Before KamLAND-2004 the bound 
on sin 2 613 from global data was dominated by the 
CHOOZ reactor experiment, together with the determi- 
nation of Am 2 ! from atmospheric data. However, in- 
cluding KamLAND-2004 the bound becomes compara- 
ble to the reactor bound. Note also that, since the reac- 
tor bound on sin 2 613 deteriorates quickly as Am 2 TM de- 
creases (see Fig. 3), the improvement is especially impor- 
tant for lower Am 2 TM values. In Fig. 3 we show the upper 
bound on sin 2 613 as a function of Am 2 TM from CHOOZ 
data alone compared to the bound from an analysis in- 
cluding solar and reactor neutrino data. One sees that, 
although for larger A;« 2 TM values the bound on sin 2 613 
is dominated by CHOOZ, for Am 2 TM < 2 x 10~ 3 eV 2 the 
solar and KamLAND data become relevant. 



FIGURE 3. Upper bound on sin 613 (1 d.o.f.) from solar and 
reactor data versus Am 2 TM . Dashed (solid) curves correspond to 
90% (3fJ) C.L. bounds, thick curves include KamLAND-2004 
data, thin ones do not. Light (dark) regions are excluded by 
CHOOZ at 90% (3<r) C.L. The horizontal line corresponds to 
the current Am 2 TM best fit value, hatched regions are excluded 
by atmospheric + K2K data at 3o~. 



Altogether, the bound on sin 2 ©13 contributes signifi- 
cantly to the overall global bound 0.047 at 3(7 for 1 d.o.f. 
As shown in [1] such an improved sin 2 613 bound fol- 
lows mainly from the strong spectral distortion found in 
the 2004 sample. 

Future long baseline reactor and accelerator neu- 
trino oscillation searches [13], as well as studies of the 
day/night effect in large water Cerenkov solar neutrino 
experiments such as UNO or Hyper-K [14] could bring 
more information on sin 2 613 [15]. With neutrino physics 
entering the precision age it is necessary to scrutize also 
the validity of the unitary approximation of the lepton 
mixing matrix in future experiments, given its theoreti- 
cal fragility [3]. 

ABSOLUTE NEUTRINO MASS SCALE 

On general grounds neutrino masses are expected to 
be Majorana [3], a fact that may explain their relative 
smallness with respect to other fermion masses. Neutrino 
oscillation data are insensitive to the absolute scale of 
neutrino masses and also to the fundamental issue of 
whether neutrinos are Dirac or Majorana particles [10, 
11]. Hence the importance of neutrinoless double beta 
decay [16]. The significance of the j3j3ov decay is given 



by the fact that, in a gauge theory, irrespective of the 
mechanism that induces j3 j3o v , it is bound to also yield 
a Majorana neutrino mass [17], as illustrated in Fig. 4. 
Quantitative implications of the "black-box" argument 
are model-dependent, but the theorem itself holds in any 
"natural" gauge theory. 

Now that oscillations are experimentally confirmed we 
know that j3/3o v must be induced by the exchange of 
light Majorana neutrinos. The corresponding amplitude 
is sensitive both to the absolute scale of neutrino mass 
as well as the two Majorana CP phases that characterize 
the minimal 3-neutrino mixing matrix [3]. Fig. 5 shows 
the estimated average mass parameter characterizing the 
neutrino exchange contribution to j3 j3o v versus the light- 
est neutrino mass. The calculation takes into account the 
current neutrino oscillation parameters from [1] and the 
nuclear matrix elements of [18] and compares with ex- 
perimental sensitivities. The upper (lower) panel corre- 
sponds to the cases of normal (inverted) neutrino mass 
spectra. In these plots the "diagonals" correspond to the 
case of quasi-degenerate neutrinos [19, 20], which give 
the largest j3/3ov amplitude. In contrast to the normal hi- 
erarchy, where a destructive interference of neutrino am- 
plitudes is possible, the inverted neutrino mass hierarchy 
implies a "lower" bound for the j3/3ov amplitude. An ex- 
ception to the rule that there is no lower bound on /3 /3ov 
in normal hierarchy models is provided by the model in 
[21]. 

Future experiments [22] will extend the sensitivity and 
provide an independent check of the Heidelberg-Moscow 
claim [23]. More information on the absolute scale of 
neutrino mass will also come from future beta decays 
searches [24] and cosmology [25]. 






FIGURE 5. j8 /3()v amplitude and current oscillation data. 



FIGURE 4. ]3j3o v and Majorana mass are equivalent [17]. 



FIGURE 6. The seesaw mechanism. 

THE ORIGIN OF NEUTRINO MASS 

Neutrino mass arise from the dimension-five operator 
U§§ where the 51/(2) ® t/(l) Higgs doublet and I is 
a lepton doublet [26]. Nothing is known from first prin- 
ciples about the mechanism that induces this operator, its 
associated mass scale or flavour structure. The neutrino 
masses that result from it once the electroweak symmetry 
breaks down are expected to be Majorana. This may ex- 
plain why neutrino masses are much smaller than those 
of the other fermions. This may happen either because 
the operator is suppressed by a large scale in the denom- 
inator, or else suppressed by a small scale in the numer- 
ator. Both ways are viable and can be made natural. 

The most popular case is the seesaw mechanism which 
induces small neutrino masses from the exchange of 
heavy states that may come from unification. Small neu- 
trino masses are induced either by heavy SU(2) (g) t/(l) 
singlet "right-handed" neutrino exchange (type I) or 
heavy scalar bosons exchange (type II), in a nomencla- 
ture opposite from the original one in [3]. The effective 
triplet seesaw term has a flavor structure different from 



the type-I term, contributing to the lack of predictivity 
of seesaw schemes. An attempt to recover predicitvity 
within the seesaw approach by appealling to extra sym- 
metries [27] is given in [20]. The model predicts max- 
imal 0?3, 013 = 0, and naturally large On, though un- 
predicted. Moreover, if CP is violated 613 becomes arbi- 
trary but the Dirac CP violation phase is maximal [28]. 
The model gives a lower bound on the absolute neutrino 
massd m v > 0.3 eV, requires a light slepton below 200 
GeV, and gives large rates for flavour violating processes. 

Amongst "bottom-up" models we mention those 
where neutrino masses are given as radiative correc- 
tions [29, 30] and models where low energy supersym- 
metry is the origin of neutrino mass [31]. The latter are 
based on the idea that R parity spontaneously break [32], 
leading to a very simple effective bilinear R parity 
violation model [33]. In this case the neutrino spectrum 
is typically hierarchical, with the atmospheric scale 
generated at the tree level and the solar scale radiatively 
"calculable" [34]. For the parameters that reproduce the 
masses indicated by current data, typically the loghtest 
supersymmetric particle decays in the detector, and its 
decay properties correlate with neutrino mixing angles, 
a test that can be made, e. g. at the LHC. 

ROBUSTNESS OF OSCILLATIONS 

The general effective model-independent description of 
the seesaw at low-energies was given in [3]. It is charac- 
terized by (n,m), n the number of SU(2) ® U(l) isodou- 
blet and m the number of SU (2) ® U{\ ) isosinglet lep- 
tons. In the mass basis a the (3,3) seesaw model has 
12 mixing angles and 12 CP phases (both Dirac and 
Majorana-type) characterizing its full 3x6 charged cur- 
rent seesaw lepton mixing matrix and non-diagonal neu- 
tral current [3]. The nontrivial structure of charged and 
neutral current weak interactions with non-unitary lep- 
ton mixing matrix is a general feature of seesaw mod- 
els [3] and lead to dimension-6 terms non-standard neu- 
trino interactions (NSI), as illustrated in Fig. 7. Such 




FIGURE 7. Flavour-changing effective operator. 

sub-weak strength eGf operators can be of two types: 
flavour-changing (FC) and non-universal (NU). In in- 
verse seesaw-type models [35] the non-unitary piece of 
the lepton mixing matrix can be sizable and hence the in- 
duced NSI may be phenomenologically important [36]. 
Sizable NSI strengths may also arise in radiative neutrino 
mass models and in supersymmetric unified models [37]. 

Non-standard physics may in principle affect neu- 
trino propagation properties and detection cross sec- 
tions [38]. In their presence, the Hamiltonian describing 
atmospheric neutrino propagation has, in addition to the 
standard oscillation part, another term //nsi 

ff NSI = ±V2G F N f (^lA- ( 2 ) 

Here +(— ) holds for neutrinos (anti-neutrinos) and e 
and e' parameterize the NSI: \f2GpNfe is the forward 
scattering amplitude for the FC process +/ — > V T +/ 
and V2GpNfe' represents the difference between +/ 
and V T + f elastic forward scattering. Here Nf is the 
number density of the fermion / along the neutrino path. 
In the 2-neutrino approximation, the determination of 
atmospheric neutrino parameters Am^ TM and sin 2 ATM 
was shown to be practically unaffected by the presence of 
NSI on down-type quarks (/ = d) [39]. Future neutrino 
factories will substantially improve this bound [40]. 

In contrast, the oscillation interpretation of solar neu- 
trino data is "fragile" in the presence of non-standard 
interactions [41], with a new "dark side" solution (with 
sin 2 9 su i ~ 0.7 [41]), essentially degenerate with the con- 
ventional one, present even after combining with data 
from reactors. On the other hand, it has been shown [42] 
that, even a small residual non-standard interaction of 
neutrinos in the e — X channel leads to a drastic loss in 
sensitivity in the 613 determination at a neutrino factory. 



It is therefore important to improve the sensitivities on 
NSI, another window of oportunity for neutrino physics 
in the precision age. 
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